Lowering Effective Coordination Promotes Adsorption of NO on Rh(lOO) and 

Rh/MgO (100) surfaces 



Raghani Pushpa^, Prasenjit Ghosh^''^, Shobhana Narasimhan^, and Stefano de Gironcoli^ 
^ Scuola Internazionale Superiore di Study Avanzati (SISSA), Via Beirut 2/4, 1-34014 Trieste, Italy 
Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur PO, Bangalore 560 064, India and 
The Abdus Salam International Centre for Theoretical Physics (ICTP), Strada Costiera 11, 1-34014 Trieste, Italy 

We have studied the adsorption of NO, and the coadsorption of N and O, on four physical and 
hypothetical systems: unstrained and strained Rh(lOO) surfaces and monolayers of Rh atoms on 
strained and unstrained MgO(lOO) surfaces. We find that as we go from Rh(lOO) to Rh/MgO(100), 
via the other two hypothetical systems, the effective coordination progressively decreases, the d- 
band narrows and its center shifts closer to the Fermi level, and the strength of adsorption and 
co-adsorption increases. Both strain and the presence of the oxide substrate contribute significantly 
to this. However, charge-transfer is found to play a negligible role, due to a cancelling out between 
donation and back-donation processes. Our results suggest that lowering effective coordination of 
Rh catalysts by strain, roughening or the use of inert substrates might improve reaction rates for 
the reduction of NO to N2. 

PACS numbers: 28.52.Fa, 68.43.Bc, 96.12.Kz, 71.20.Be 



I. INTRODUCTION 

Though catalysts are crucial to the operation of many 
industrial and commercial processes, our understanding 
of the factors that make a good catalyst is still incom- 
plete. Both the electronic and geometric structures of 
the catalyst are known to be important, and in recent 
years, it has become evident that coordination number 
can have a large effect on catalytic activity. The gen- 
eral understanding is that lower coordination leads to 
higher activity, though there arc exceptionSfi Industrial 
catalysts usually consist of small particles supported on 
a substrate. The presence of the substrate introduces 
additional factors, such as charge transfer between the 
catalyst and the support, a change in the local environ- 
ment of catalyst atoms, and geometrical strain imposed 
by the presence of the substrate. 

In this paper, we have studied some of these issues by 
performing calculations to study the adsorption of NO, 
as well as the coadsorption of N and O, on a variety of 
realistic and hypothetical Rh(lOO) surfaces. These con- 
figurations constitute the initial and final states in the 
dissociation of NO, which is a crucial step in the reduc- 
tion of nitrogen monoxide to nitrogen, for example in 
three-way catalysts in automobiles. Rh is perhaps the 
best catalyst for this process^ but it is also by far the 
most expensive precious metal. It is therefore desirable 
to gain a better understanding of what makes a good cat- 
alyst, so as to guide one in developing new catalysts that 
use less or no Rh. 

Towards this end, and in an attempt to gauge the con- 
tributions made by the chemical nature of the substrate, 
as well as strain, we have studied adsorption and coad- 
sorption on four kinds of (100) surfaces: (i) the surface of 
an unstrained Rh crystal; (ii) the surface of a Rh crystal 
that has been stretched (expanded) so as to be commen- 
surate with an MgO(lOO) surface; (iii) a monolayer of 
Rh on MgO(lOO), with both the Rh and the substrate's 



in-plane spacing fixed to be that of Rh(lOO); and (iv) a 
monolayer of Rh on MgO(lOO), with all in-plane spac- 
ings fixed to be that of MgO(lOO). Note that straining 
systems (i) and (iii) leads to systems (ii) and (iv) respec- 
tively, while changing the substrate from Rh to MgO 
converts systems (i) and (ii) to systems (iii) and (iv) re- 
spectively. Qualitatively, it seems apparent that as one 
proceeds from system (i) to system (iv), the surface Rh 
atoms become progressively less coordinated; this idea is 
put on a quantitative footing further below! . 

Another factor - magnetism - seems likely to play a 
role in these systems. Bulk Rh is "almost" ferromag- 
netic. There has been a long-running controversy over 
whether or not the Rh(lOO) surface is magneticj^i^i^ii 
however it seems clear that Rh monolayers and clusters 
are magnet ici^ i^dOd^ This raises the question of what 
possible role magnetism may play in the operation of Rh 
catalysts; this issue is also dealt with in this paper. 

Though the adsorption of NO on Rh(lOO) has been 
studied by previous authorS )^'^'^^ we are not aware of 
any systematic programme of calculations that is similar 
in spirit to ours. 



II. METHOD 

Our calculations have been performed using ah initio 
density functional theory, using the PWscf package of the 
Quantum-ESPRESSO distributioujH The spin-polarized 
version of the Kohn-Sham equations were solved using 
ultrasoft pseudopotentialsi^ and a plane-wave basis with 
a cut-off of 30 Ry. Exchange and correlation effects 
were treated using the Generalized Gradient Approxi- 
mation (GGA) in the form suggested by Perdew, Burke 
and Ernzerhofji^ In order to improve convergence, a 
Methfessel-Paxton smearingi^ with a width of 0.03 Ry 
was used. 

Most results for the Rh(lOO) and Rh/MgO(100) slabs 
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FIG. 1: Schematic diagrams showing the geometries consid- 
ered for adsorption studies. The grey, white and black circles 
represent Rh, O and N atoms respectively. The NO occupies 
different sites, and is oriented differently, in the six cases de- 
picted here. (a),(b) and (c) show top views of configurations 
where the NO sits horizontally on the Rh substrate, while (d), 
(e) and (f) show top views of configurations where the NO sits 
vertically on the substrate; in these, the N atoms are not vis- 
ible, as they sit directly below the O atoms. The picture (g) 
shows a side view of configuration (e). 




were obtained by using a (1 x 1) asymmetric supercell 
containing four Rh layers, of which the outermost two 
(towards the adsorbate) were allowed to relax, while the 
inner two were kept fixed at the appropriate bulk separa- 
tion. However, some tests were also performed with sym- 
metric slabs containing eight layers of Rh atoms. Bril- 
louin zone integrations for such (1x1) surface cells were 
carried out using a (12x12x1) Monkhorst-Pack mesh.— 

The adsorption and co-adsorption studies were carried 
out using slabs with larger unit cells, [viz., (2x2) and (2x 
3)], together with corresponding k-point meshes [(6 x 6 x 
1) and (6x4x1) respectively]. For systems (i) and (ii), we 
found that a slab containing four Rh layers was sufficient 
to give well-converged adsorption energies; however, for 
one particular adsorption geometry, we found it necessary 
to use a slab with five Rh layers. All the results presented 
below for adsorption and co-adsortpion on systems (iii) 
and (iv) were obtained using a slab containing one layer 
of Rh atoms over four MgO layers. 

In adsorption studies, six different geometries were 
considered; these are shown in Fig. [U while Fig. [2] shows 
the three different geometries for the co-adsorption stud- 
ies. These two figures show (schematically) the geome- 
try prior to relaxation; after relaxation (in the absence of 
symmetry constraints) the geometries remained roughly 
similar, though bond lengths changed. 

As we will show in the next section, in cases (ii), (iii) 
and (iv), the Rh(lOO) surface was found to be magnetic. 
In such cases, one needs to worry about the possibil- 
ity that the bottom surface acquires a magnetic moment 
and is thus no longer representative of a bulk-like layer. 
To mitigate this, in these cases we deposited a layer of 
H atoms on the bottom surface, which has the effect of 
quenching magnetization on that surface. 



FIG. 2: Top view of the geometries considered for N and O 
adsorption on Rh(lOO) surfaces. The grey, white and black 
circles represent Rh, O and N atoms respectively. In config- 
uration (a), N and O are at nearest neighbour hollow sites, 
in (b) N and O are at diagonally opposite hollow sites. In 
configuration (c), N and O are at the bridge sites between 
two Rh atoms. The white rectangle indicates the 2x3 unit 
cell used for co-adsorption studies. 



(a) (b) (c) 




III. RESULTS 



A. Preliminary Tests: bulk Rh, bulk MgO and 
gas-phase NO 

We obtained a value of 3.85 A for the lattice constant of 
bulk Rh, which is in excellent agreement with the exper- 
imental value of 3.80 and previous theoretical values 
of 3.87 Af^^ Our value implies that the distance between 
nearest-neighbor (NN) Rh atoms on the (100) surface is 
2.73 A, while the interlaycr distance between bulk-like 
layers in a Rh(lOO) slab is 1.93 A. 

For MgO, we obtained a lattice constant of 4.25 A, 
which is identical to the value obtained in previous 
calculations,— and in good agreement with the exper- 
imental value of 4.21 Ai2i Note that this would imply 
that if a Rh monolayer were to be deposited commensu- 
rately with an MgO(lOO) substrate, the NN Rh-Rh dis- 
tance within the monolayer would be increased to 3.00 
A, which corresponds to a strain of 9.9%. 

For NO in the gas phase, we obtained a binding en- 
ergy of 7.13 eV, and an N-O bond length of 1.17 A . For 
comparison, the experimental values are 6.5 eV and 1.15 
A respectively!^ 

In all these cases, it can be seen that our results are 
in good agreement with both experiments and previous 
calculations, lending support to the validity of our ap- 
proach. 



B. Case (i) - Unstrained Rh(lOO) surface 

For the surface energy of the clean and unstrained 
Rh(lOO) surface [case (i)], we obtain a value of 1.12 eV 
per surface atom, which agrees exactly with the value 
obtained in a previous calculatiouji^ and is also in rea- 
sonably good agreement with the experimental value of 
1.27 cV per surface atom<^ We find that the first in- 
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terlayer spacing di2 is contracted by 3.5% with respect 
to the bulk interlaycr spacing of 1.93 A; this is similar 
to the contraction of 4.0% found in an earlier study,— 
but more than the contraction of 1.4% ± 1.4% reported 
experimentally)^ We find that the next two interlayer 
spacings, c?23, and ^34, are both expanded by 0.77%; the 
next interlayer spacing ^45 is very close to the bulk in- 
terlayer spacing. We find that the surface is not ferro- 
magnetic, which is in agreement with several previous 
studies,— though it disagrees with some reports in the 
literature^i^ 

Following previous authors ^21 we define an effective co- 
ordination number for the surface as: 

= (1) 

where, p'jfjiiR) is the atomic charge density of an isolated 
Rh atom as a function of R, the distance from the nu- 
cleus, and the sum is taken over all the nearest neighbour 
atoms j around a surface Rh atom i. Rij is the distance 
between atoms i and j, and Rbuik is the NN distance 
for bulk Rh. For case (i), we find = 8.49. This is 
slightly increased with respect to the nominal surface co- 
ordination (the number of nearest neighbor atoms for a 
surface atom) of 8, due primarily to the contraction of 
di2 relative to the bulk interlayer spacing. 

We define the adsorption energy as Eads = 
ENO:B.h(WQ) - ERh(wQ) - E^j^q, where E is the value ob- 
tained from ah initio calculations for the total energy 
of the corresponding configuration. In agreement with 
previous results ^iii^ we find that the most favorable ad- 
sorption geometry on an unstrained Rh(lOO) slab is the 
"vertical bridge" (VB) [corresponding to FigsllUe) and 
(g); the adsorption geometry is also shown in Fig. [3Ua)]. 
With a (2 X 2) unit cell (i.e., a coverage of 1/4 ML), we 
obtain Eads = -2.59 eV in this geometry; and the N-0 
bond length is increased slightly from the gas-phase value 
of 1.17 A to 1.20 A. The next most favorable adsorption 
geometry is the "horizontal hollow" (HH) [see FiglT](a)] 
with Eads = -2.47 eV, and an N-0 bond length of 1.31 
A in the latter geometry. The details of the adsorption 
geometries for both VB and HH are given in Table H] All 
these values are in excellent agreement with a previous 
studyfi^ 

We have also investigated the geometry and energet- 
ics of the co-adsorption of N and O on this surface. We 
find that both N and O prefer to occupy four-fold hollow 
sites, as has also been reported in a previous paper In 
Table Ull we show how the co- adsorption energy, defined 

as Ecoads = Ei^+0:Rh(wo) " £'fl/i(ioo) " ^wo' varics with 
adsorbate coverage and adsorption geometry. We find 
that the N and O atoms prefer to occupy next-nearest- 
neighbor hollow sites over nearest-neighbor hollow sites. 
Also, the magnitude of Ecoads increases as the coverage 
is decreased. These observations suggest that there is a 
repulsive interaction between the N and O atoms after 
dissociation of NO. We note that in a previous study a 
coadsorption energy of -3.50 eV was obtained at a cover- 




(a) (b) 




(c) (d) 

FIG. 3: (In color online). The most stable configurations for 
adsorbed NO in (a) case (i), (b) case (ii), (c) case (iii), and (d) 
case (iv). In cases (i) and (ii), vertical configuration of NO 
at a bridge site is the most stable; whereas in cases (iii) and 
(iv), horizontal configuration at the hollow site is the most 
stable. Small blue and red, and large grey and turquoise 
spheres represent N, O, Rh and Mg atoms respectively. 



age of 1/16 ML, which is in keeping with our conclusions 
here. 



C. Case (ii) - Stretched Rh(lOO) surface 

When an in-plane expansion of 9.9% is imposed on 
a (1x1) Rh(lOO) slab, we find that the bulk interlayer 
spacing (w ^45) is reduced by 11.3% with respect to the 
unstrained case and becomes equal to 1.71 A. Upon al- 
lowing all interlayer distances to relax in an eight-layer 
slab, the first three interlayer spacings (ii2, ^23 and ^34 
are reduced to 1.59, 1.76 and 1.73 A respectively. These 
are contracted significantly, with respect to the bulk in- 
terlayer distance in unstrained Rh(lOO), by 17.3, 8.9 and 
10.5 % respectively. The net effect of having longer in- 
tralayer distances but shorter interlayer distances at the 
surface is that now the effective coordination Uf. = 6.65, 
i.e., as a result of stretching, the effective coordination 
has decreased on going from case (i) to case (ii). 

In Fig. m we show the layer-resolved magnetization 
per atom of a symmetric (1x1) stretched eight-layer 
slab (stars). Interestingly, the stretched Rh(lOO) sur- 
face is now magnetic with surface magnetization of 1.0 
/is /atom. In order to reduce computational load, it is 
a common practice to consider asymmetric slabs where 
some of the bottom layers are fixed at their bulk posi- 
tions. The circles in Fig. 5] show the magnetization per 
atom in a (1x1) four-layer slab with the bottom two lay- 
ers fixed at 1.71 A interlayer spacing and the top two 
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TABLE I: The Rh-Rh bond length (dRh-Rh), the Rh-N bond length (dflh-jv), the Rh-O bond length {dnh-o) and the N-0 
bond length {(In-o) of NO on different adsorption sites for the four cases. The number given below each case denotes the 
Rh-Rh in plane distance for the clean surfaces. RHm and Rho indicate the distances between two Rh atoms bonded to N and 
O atoms respectively. All distances are given in A. 



System 


Site and geometry 


dRh-Rh 


dRh-N 


dRh-o 


d]v-o 


case (i) 


VB 


2.74 


1.96 


- 


1.20 


(2.73) 


HH 


2.75 (Rhiv) 
2.79 (Rho) 


1.98 


2.21 


1.31 


case(ii) 


VB 


2.70 


1.96 




1.20 


(3.00) 


HH 


3.13 (Rhiv) 
3.21 (Rho) 


1.98 


2.22 


1.32 


case (iii) 


VB 


2.75 


1.93 




1.20 


(2.73) 


HH 


2.70 (Rhiv) 
2.95 (Rho) 


1.95 


2.15 


1.33 


case(iv) 


VB 


2.56 


1.94 




1.21 


(3.00) 


HH 


2.66 (Rhjv) 
2.59 (Rho) 


1.95 


2.18 


1.34 



TABLE II: Dependence of coadsorption energy Ecoads on cov- 
erage and adsorption geometry. The first row corresponds to 
a coadsorption of N and O in a 2 x 2 unit cell in next-nearest- 
neighbor hollow sites, while the second and third rows (a) 
correspond to results for the geometries depicted in Figs[2la) 
and (b) respectively. Here, C = coverage of N = coverage of 
O. 



Coverage C 


Ecoads (eV) 


rjv-flh(A) 


ro-Rh (A) 


1/4 ML 
1/6 ML (a) 
1/6 ML (b) 


-3.20 
-3.10 
-3.34 


2.03 
1.99, 2.10 
2.03, 2.04 


2.14 
2.04, 2.35 
2.14, 2.17 



layers relaxed. Wc find that only in the top two lay- 
ers is the magnetization close to the one in the layers of 
the symmetric slab, and a "spurious" magnetization as 
large as 1.2 yu^/atom is present on the frozen-geometry 
bottom surface. In order to quench this spurious magne- 
tization, wc considered also an asymmetric four-layer Rh 
slab where H atoms were adsorbed on the bottom sur- 
face. The results are again reported in Fig. [5] as squares, 
and wc now find that the magnetic moments on all layers 
are very similar for the four-layer and eight-layer slabs, 
giving us confidence that adsorption energies obtained 
with the former will be accurate. 

The interlayer distances are also very similar in the 
eight-layer symmetric slab and in the four-layer slab with 
H adsorbed on the bottom surface, but are (slightly but 
perceptibly) different on the four-layer slab without H on 
the bottom surface. 



Adsorbing H on the bottom surface also turns out to 
be crucial in correctly predicting the stable adsorption 
geometry for NO. In the absence of H, it is found that 
the VB and HH configurations are degenerate, with both 
having an adsorption energy Eads = -2.83 eV; however, 
when H is adsorbed on the bottom surface, the VB (with 
Eads ~ -2.86 eV) is found to be clearly favored over the 
HH (with Eads = -2.71 eV); for the reasons stated above, 
we believe that the set of numbers obtained with H atoms 
on the bottom surface is more to be trusted. This stable 
adsorption geometry is depicted in Fig. [3jb). 

In order to gauge the effects of magnetization on ad- 
sorption energies, we also performed non-spin-polarized 
(NSP) calculations. Wc found that the magnitude of 
Eads increases in the absence of magnetization; how- 
ever, the difference between the VB and HH geometries 
is maintained: with NSP calculations, the former gives 
Eads = —2.96 eV, while the latter gives Eads = —2.82 
eV. 

Wc found that NO adsorption leads to a significant dis- 
tortion in the position of Rh atoms (see Table [l|. When 
NO is adsorbed in the VB configuration, the Rh-Rh dis- 
tance for the two Rh atoms bonded to N is reduced sig- 
nificantly, from 3.00 A to 2.70 A. However, for the HH 
configuration, the distances between the two Rh atoms 
bound to N and O increase to 3.13 and 3.21 A respec- 
tively. The internal bond length in NO is increased in 
all cases; to 1.20 A for the VB and to 1.32 A for the HH 
configuration. Note that this also suggests that it might 
be easier to break the N-0 bond in the HH configuration. 

Next, we co-adsorb N and O on this surface. The co- 
adsorption geometries considered by us are shown in Fig. 
[51 When we start from a configuration where N and O are 
at NN hollow sites (see Figl^Ja)), the O atom moves away 
from N, and sits at the bridge site which is equidistant 




Layer number 



FIG. 4: Magnetization of the stretched Rh(lOO) surface [case 
(ii)] for an eight- layer symmetric slab (stars), a four-layer 
asymmetric slab (cirlces) and a four-layer asymmetric slab 
with H adsorbed at the bottom of the slab (squares). 



from the N atom and its image. This configuration rep- 
resents the most favorable co-adsorption geometry, with 
Ecoads = —3.47 eV. The second most favourable config- 
uration is one where both N and O sit at NNN hollow 
sites (Fig. HJ^b))); here Ecoads = —3.33 eV. Unhke in 
case (i), the stretched Rh(lOO) surface also has a local 
minimum when both the atoms sit at bridge sites (see 
Fig. (He)). This is the configuration that the system as- 
sumes just after the dissociation of NO, and corresponds 
to Eads = —3.24 cV. As in case (i), the magnitude of 
Eads is increased by ^ 0.3 eV in the first two adsorption 
geometries upon performing an NSP calculation; how- 
ever, in the third case, Eads is essentially unchanged. It 
is also worth noting that changes in interatomic distances 
are slight when comparing spin-polarized (SP) and NSP 
calculations. 



D. Case (iii) - Monolayer of Rh on MgO(lOO), at 
Rh(lOO) lattice constant 

The third case we consider is again an artificial one. 
We consider first an MgO(lOO) substrate that has been 
strained (compressed) in-plane, so that it has the same 
in-plane lattice constant as in case (i). As a result of 
this in-plane compression, we find that the first two in- 
terlayer spacings are both increased by 8.15% relative to 
the interlayer spacing in unstrained MgO. There is also a 
noticeable rumpling of Mg-0 layers near the surface. We 
define the rumpling as Az = {zo—ZMg)/do, where zq and 
ZMg are the z-coordinates (normal to the surface plane) 
of O and Mg respectively, and do is the bulk interlayer 
spacing. Then, the outermost MgO layer has Az ~ 2.0%, 
while the second MgO layer has Az = —0.1%. Note that 
the rumpling is opposite in the two outermost layers: in 
the topmost layer, oxygen atoms are displaced further 
away from the substrate than Mg atoms, whereas for the 
second layer, the reverse is true. 

When a monolayer of Rh is deposited pseudomor- 
phically on this compressed MgO substrate, it binds 
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FIG. 5: Charge transfer on placing a monolayer of Rh atoms 
on an MgO substrate for (a) case (iii) and (b) case (iv). Red 
and dark blue lobes represent regions where the charge has in- 
creased and decreased respectively. Red and turquoise spheres 
represent O and Mg atoms respectively, the grey Rh atoms 
are not visible because they are surrounded by charge-transfer 
lobes. Note that both the Rh overlayer and the O atoms in 
the topmost Mg-0 layer display both red and blue lobes. 



with a binding energy Ehin of -3.59 eV per Rh atom. 
Here, the binding energy has been defined by Ei,in = 
Em.MgO - Emqo - E^^^, where Enh-.MgO is the total 
energy of the Rh/MgO slab, EMgO is the energy of the 
MgO slab alone, and E^^^ is the energy of an isolated 
Rh atom in the gas phase. The Rh atoms sit atop the 
oxygen atoms of the outermost MgO layer, with a Rh-0 
separation of 2.22 A. The deposition of the monolayer 
changes both the rumpling and the relaxation of the out- 
ermost MgO layers. The first two layers now display a 
rumpling of -4.2 % and 0.7% respectively; note that the 
deposition of Rh has actually reversed the direction of 
rumpling. There is also a change in interlayer spacings, 
with the distance between the two outermost MgO layers 
now expanded by 9.55% relative to the interlayer distance 
in unstrained MgO, while the expansion of the next in- 
terlayer spacing is now reduced slightly, from the value of 
8.15% (in the absence of the Rh adlayer) to 8.05%. The 
Rh monolayer in this case is found to be magnetic, with 
a magnetic moment of 1.5 /.is /atom. 

In principle, placing metal atoms on an oxide substrate 
can lead to charge transfer. In Fig. [5]Ja), we have plot- 
ted the redistribution of charge upon placing a layer of 
Rh atoms on an MgO substrate. In this figure, red and 
dark blue regions represent areas where the charge has 
increased and decreased respectively. It is interesting to 
note that one observes both red and dark blue lobes on 
both the Rh layer and the topmost Mg-0 layer; this is be- 
cause of the simultaneous presence of donation and back- 
donation between the adsorbate (Rh) and the substrate 
(MgO), as is well-known for such systems. Interestingly, 
these two processes essentially cancel out for this sys- 
tem, the net charge transfer, as calculated by Lowdin 
population analysis, is only 0.015 electrons per Rh atom 
between the Rh overlayer and the MgO substrate. 

The effective coordination is now lowered to 5.51. 
Note that Up is lower for both cases (ii) and (iii) relative 
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to case (i); however, the lowering in case (iii) (which is 
due to the replacement of the Rh substrate by an MgO 
support), is more significant than in case (ii) (where it 
is due to stretching the system to the lattice constant of 
MgO). Similarly, the magnetic moment is larger in case 
(iii) than in case (ii). 

Upon adsorbing NO on this compressed Rh/MgO(100) 
system, we find that the most favorable adsorption geom- 
etry is the HH [depicted in Fig.[3fc)], with Eads ~ —3.38 
eV, which is larger in magnitude than the adsorption en- 
ergy of -3.17 eV for the VB configuration. Note that this 
is in contrast to the most favored adsorption geometry 
for both cases (i) and (ii); it is now favorable for both N 
and O to bind to Rh atoms. 

For the VB case, the distance between the two Rh 
atoms bound to the N (RIiat) increases slightly from 2.73 
A to 2.75 A. On the contrary, for the HH geometry, the 
distance between the same pair of Rh atoms decreases to 
2.70 A, while that between the Rh atoms attached to the 
O atom (Rho) increases to 2.95 A; this is because, as a 
result of strong binding between NO and the Rh atoms, 
the Rh^v pair is pulled out of the surface. 

Moreover, upon adsorbing NO, the magnetization of 
the Rh atoms is quenched. For the HH NO adsorption 
geometry, this reduction is as follows: for the two Rh 
atoms to which N is attached, the reduction is by 74%, 
for the two Rh atoms to which O is attached, it is by 
45%, and for the two Rh atoms which are furthest from 
NO, the reduction is only by 9%. 

Once again, we find significantly enhanced binding 
upon performing an NSP calculation, with Eads = —4.07 
eV and -3.82 eV for the HH and VB geometries respec- 
tively. 

The co-adsorption energies are -4.64 -3.32 eV when 
the N and O atoms are in the NN hollow (Fig.[2Ka)) and 
bridge (Fig. [IJc)) sites respectively. 



E. Case (iv) - Monolayer of Rh on MgO(lOO), at 
MgO lattice constant 

Upon putting a pseudomorphic layer of Rh atoms on an 
unstrained MgO (100) substrate, we find that Rh atoms 
again preferentially occupy the sites atop O atoms. The 
binding energy for this configuration, defined as in case 
(in), is found to be Etin = —3.27 eV per Rh atom; this 
is somewhat less in magnitude than the value of -3.99 
eV obtained in a previous studyi^ The Rh-0 distance 
is found to be 2.25 A, which is larger than the value of 
2.10 A reported earlicriiSll Both the Mg and O atoms in 
the topmost layer move outwards, away from the sub- 
strate, so that the mean interlayer distance between the 
top two MgO layers is increased by 0.73% relative to the 
bulk interlayer spacing. However, there is a considerable 
rumpling (Mg atoms are higher up on the surface than 
oxygen atoms), with Az = —4.0%, of the topmost MgO 
layer; this is in good agreement with one earlier reported 
value of -4.3%iiS but somewhat less than the value of 



-6.1% obtained in another studyi^S The rumpling is re- 
duced to +1.0% in the MgO layer below this one. Note 
that for both layers the rumpling is reversed with respect 
to that observed for a bare MgO(lOO) surface. 

For a Rh atom in case (iv), we obtain = 3.59. Note 
that this is significantly lowered with respect to the value 
of 8.49 obtained in case (i), due to a combination of two 
factors: the replacement of substrate Rh atoms by Mg 
and O atoms, and the stretching to the MgO lattice con- 
stant. 

Once again, we find that the surface Rh atoms are 
magnetic, with a moment of 1.5 ^s/atom. 

In Fig. [51^b), we present a plot of charge transfer for 
the Rh/MgO system. From this figure, it is evident that, 
as in case (iii), there is both donation and back-donation 
between the overlayer Rh atoms and the O atoms in the 
MgO substrate. Once again, these two processes effec- 
tively cancel out, and there is a net charge transfer of 
only 0.018 electron per Rh atom from the oxygen to the 
Rh atom. 

Next, we study the adsorption of NO on this 
Rh/MgO(100) system. As in case (iii) above, we find 
that it is now most favorable for the NO atom to lie hor- 
izontally on the surface: Eads = —4.08 eV and -3.87 eV 
for the HH and VB geometries respectively. In Fig. [3] (d) 
we have depicted the lowest-energy adsorption geometry, 
corresponding to the HH configuration. 

With the adsorption of NO on 1 ML of Rh on MgO, the 
Rh atoms to which the NO is attached come closer to- 
gether for both the HH and VB geometries (see Table 
We speculate that this may be because Rh-Rh bonds are 
under considerable tensile stress when deposited pseudo- 
morphically on MgO; the adsorption of NO breaks the 
symmetry and allows the Rh atoms bonded to NO to 
come closer together. 

As in case (iii), upon adsorbing NO, the magnetiza- 
tion of the Rh atoms is quenched. Out of the six Rh 
atoms present in the unit cell, the spin polarization of 
each of the two Rh atoms to which N is attached is neg- 
ligible, while those to which O is attached is reduced by 
about 81% and the spin polarization of the two Rh atoms 
which are furthest from the NO is reduced by about 28%. 
Similar effects have also been observed previously, by us 
for NO adsorbed on small Rh clusterSf^^ and by Hass et 
al. in their studies of NO adsorption on a hypothetical 
monolayer of Rh atoms Apart from the fact that the 
latter group of authors worked with a monolayer of Rh 
(i.e., there was no substrate), there are other differences 
between our calculations and theirs: they did the calcu- 
lations for NO adsorbing in the VB geometry alone, fixed 
the Rh-Rh distance at the value for bulk Rh, and did not 
relax the coordinates of NO. While they found that the 
magnetic moments on all the atoms of the monolayer are 
negligible, we find that the magnetism of those Rh atoms 
which are attached to the NO molecule is very strongly 
suppressed, while in the other Rh atoms there is a slight 
reduction. 

As in the previous cases, the magnitude of Eads in- 
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creases upon performing an NSP calculation; one obtains 
values of -4.38 and -3.93 eV for the HH and VB configu- 
rations respectively. 

Upon co-adsorbing NO on this Rh/MgO slab, wc find 
that it is preferable for both the N and O to be at NN 
hollow sites rather than at NN bridge sites; the former 
leads to a co-adsorption energy of -4.74 eV, while the 
latter leads to a co-adsorption energy of -4.10 eV. 



IV. DISCUSSION OF TRENDS 

We have seen that as we go from case (i) to case (iv), 
the effective coordination Ue decreases progressively and 
significantly. This decrease in effective coordination is 
accompanied by (and, in our interpretation, causes) a sig- 
nificant increase in the strength of adsorption of NO and 
co-adsorption of N and O. Our main results are encapsu- 
lated in Fig.[6l In the panel on the left-hand-side, we have 
presented our results for the adsorption energy of NO on 
these systems. The open circles represent results for the 
vertical bridge geometry, and the filled circles those for 
the horizontal hollow geometry. The following features 
are evident from this graph: (a) for a given geometry, 
the magnitude of Eads increases monotonically as rig is 
decreased; (b) the slope of this graph is larger for the HH 
than for the VB; this can be rationalized as being due to 
the fact that for the HH, both the N and O are bonded 
to Rh atoms, whereas for the VB, only N is bonded to 
Rh atoms; (c) as a consequence of this the HH geome- 
try becomes more favorable at lower Ue', (d) the slope of 
the lines connecting cases (i) and (ii) is less than that 
connecting (ii), (iii) and (iv). This last observation can 
be attributed to the fact that in case (i) the substrate is 
non-magnetic, whereas in case (ii) it is magnetic. In the 
absence of magnetism, adsorption is stronger, and the 
points corresponding to cases (ii), (iii) and (iv) would 
be shifted downwards, and one would have obtained a 
roughly straight line connecting the points from the four 
cases considered by us; since these three points have, 
however, been shifted upwards by the presence of mag- 
netism, the slope connecting (i) and (ii) is reduced. The 
main conclusion from this graph is that the adsorption 
energy scales more-or-less linearly with effective coordi- 
nation; we note further that the Rh(lll) surface has a 
larger (9.15) and weaker Eads (-2.18 cV) than all four 
cases considered by us. 

Similar conclusions about co-adsorption can be drawn 
from the right-hand-side panel, where the co-adsorption 
energy is found to vary monotonically with rig . The vari- 
ation is less linear for this case, presumably because of 
changes in co-adsorption geometry. 

The adsorption energy decreases as we go from (i) to 
(iv) due to two reasons: (a) the strength of adsorption in- 
creases as the coordination is decreased, and (b) the hor- 
izontal hollow geometry, which leads to stronger binding 
to the substrate and a weaker and longer NO bond, be- 
comes more favored at low effective coordination. How- 
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FIG. 6: Variation of (a) adsorption energy of NO, Eads, and 
(b) co-adsorption energy of N and O, with effective coordi- 
nation Tie. The open and filled circles in (a) represent the 
vertical bridge and horizontal hollow adsorption geometries 
respectively. The squares in (b) represent the co-adsorption 
energy for the lowest energy configuration of all geometries 
considered by us. 



ever, these effects would have been even more marked 
if the substrate were to remain non-magnetic: in every 
case, we have seen that NSP calculations (where mag- 
netism is suppressed) point to stronger adsorption than 
SP calculations. 

For the VB geometry, the N-0 bond length is found 
to be 1.2 A, whereas for the HH geometry, the bond 
length is increased to ~ 1.3 A, suggesting that it should 
be easier to break the N-0 bond in the latter case. 

Our results show that strain and the presence of the 
oxide substrate contribute to the increased strength of 
adsorption primarily through their effect on rie, and that 
both effects contribute to roughly the same extent. In the 
case of the magnesia substrate considered by us, charge 
transfer plays a negligible role, since the donation and 
back-donation mechanisms essentially cancel out. This 
may not be true for more active oxide substrates, such 
as titania and ceria. 

It seems intuitively obvious that when Rh atoms have 
a lower effective coordination, they will bind adsorbates 
more strongly, thus making it easier to break bonds 
within the adsorbate. An alternative and equivalent way 
of looking at this effect is to consider the effect of lower 
coordination on the density of states. In Fig. [7l we show 
how the spin-polarized d-band density of states of the 
surface Rh atoms changes as we go from case (i) to case 
(iv) - it can be seen that the lowering of results in a 
progressive narrowing of the rf-band. The d-bandwidth 
and position of the d-band center have been shown to be 
a good predictor of catalytic activity— 

In Fig. [5]Ja) we show how the position of the d-band 
center shifts with effective coordination. We see that as 
we go from case (i) to case (iv), the reduction in rie is 
accompanied by an approximately linearly proportional 
shift in the d-band center, bringing it closer to the Fermi 
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FIG. 7: (Color online) Rh d-band projected density of states 
(DOS) for the four cases. Positive and negative values denote 
the DOS for spin-up and spin-down electrons respectively. 



level. In Fig [Hb) we show that the adsorption energy 
varies in a monotonic (and approximately linear) way 
with the d-band center. 



V. SUMMARY AND CONCLUSIONS 

We have studied the adsorption of NO and the co- 
adsorption of N and O on strained and unstrained Rh 
surfaces with and without the presence of an MgO sub- 
strate. Both strain and placing a monolayer of Rh atoms 
on the oxide substrate lead to a significant lowering in the 
effective coordination of surface Rh atoms; doing both 
[i.e., placing a monolayer of Rh atoms pseudomorphically 
on an MgO(lOO) substrate] leads to the largest decrease 
in effective coordination rig. Further, both strain and 
the presence of the substrate (either separately or to- 
gether) have the effect of making the surface Rh atoms 
magnetic. Every decrease in Ue is accompanied by (and, 
presumably, causes) a decrease in the d-bandwidth of Rh 
atoms, a shift of the d-band center towards the Fermi 
level, and a strengthening of the adsorption of NO and 
the CO- adsorption of N and O. 

We note that the effective coordination is a quan- 
tity that can be very simply computed, especially if the 
structure is known - it is not even necessary to perform 
an ab initio density functional theory calculation in or- 
der to compute it. Therefore it can serve as a simple 
guide or rule-of-thumb in order to design systems where 
the strength of adsorption or co-adsorption takes on a 
desired value. 

Thus, lowering the effective coordination seems to be 



a good strategy to increase the strength of adsorption 
and co-adsorption, and thus, conceivably, lower the bar- 
rier to dissociation of the NO bond. One can think of 
several ways of reducing effective coordination: e.g., by 
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FIG. 8: Variation of (a) the position of the d-band center 
{Ed) with effective coordination, and (b) the adsorption en- 
ergy Eads with Ed- In both cases, a monotonic and approx- 
imately linear relationship is observed. The points plotted 
correspond to the lowest-energy configuration for each case. 



using rough surfaces, by placing Rh atoms on an inert 
oxide substrate, and by using Rh nanocatalysts where 
the Rh particles are sufficiently small so as to be signifi- 
cantly under-coordinated. However, most of these strate- 
gies to reduce Ug also favor magnetism, which competes 
with bondingi^ This suggests that an optimal value of 
Ue should be aimed for, where the coordination is low 
enough so as to favor a significant strengthening of Rh- 
NO bonds and weakening of the N-0 bond, yet the ad- 
verse effects of magnetism are not too evident. 

We are in the process of computing dissociation barri- 
ers, in order to verify our expectation that lowering ef- 
fective coordination will also weaken the N-0 bond and 
thus catalyze the reduction of NO to N2. 
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